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We determine the location A c of the mobility edge in the spectrum of the hermitian Wilson 
operator on quenched ensembles. We confirm a theoretical picture of localization proposed for 
the Aoki phase diagram. When A c > we also determine some key properties of the localized 
eigenmodes with eigenvalues |A| < A c . Our results lead to simple tests for the validity of simulations 
with overlap and domain-wall fermions. 

PACS numbers: 11.15.Ha, 12.38.Gc, 72.15.Rn 



Localization of electronic wave functions is a familiar 
phenomenon in disordered systems . Recently we con- 
jectured that a similar phenomenon takes place in lat- 
tice QCD with Wilson fermions: In an ensemble of gauge 
configurations, the low-lying eigenmodes of the hermitian 
Wilson operator can be localized, up to some mobility 
edge above which they become extended. This obser- 
vation has two important applications. First, it helps 
resolve a paradox in the quenched theory with negative 
bare mass mo, where simulations with two valence quarks 
have discovered a condensate that breaks the isospin 
symmetry in regions without Goldstone bosons Q . Sec- 
ond, there are important implications for large-scale sim- 
ulations of QCD with domain-wall |j] and overlap |5j 
fermions. Both of these formulations are based on Wil- 
son fermions with negative mo. Quenched as well as un- 
quenched calculations with these fermions will thus be 
sensitive to the spectrum of the Wilson operator. It turns 
out that an understanding of the localization properties 
is important for ensuring chirality and locality. 

In two-flavor QCD with Wilson fermions, part of the 
"supercritical" region (—8 < mo < 0) is the so-called 
Aoki phase 0, where a pion condensate (^3) breaks both 
parity and isospin symmetry. Inside the Aoki phase one 
pion is massive, whereas the other two pions are Gold- 
stone bosons. Outside the Aoki phase (e.g., for weak 
coupling, away from the critical values mo = 0, —2, . . .) 
all pions are massive. 

It is this supercritical massive phase that presents 
the conundrum in the quenched theory. The word 
"quenched" here can refer to any ensemble of gauge con- 
figurations generated without the Wilson-fermion deter- 
minant, whether a pure gauge ensemble or an ensem- 
ble with dynamical fermions of some other type, such as 
domain- wall fermions jl4| . Both the massless and the 
massive supercritical phases support a non-zero density 
of near-zero eigenmodes, and hence a condensate through 
the Banks-Casher relation (^3} = 27rp(0) 7], where p(X) 
is the spectral density of the hermitian Wilson operator 
H\y . Why are there no Goldstone bosons? In Ref. |2J we 



proposed a detailed physical picture which resolves this 
puzzle. A key observation (first made in Ref. 8) is that, in 
a quenched system, localization provides an alternative 
to the Goldstone theorem. In this Letter we present nu- 
merical evidence supporting and illustrating this picture. 
The implications for domain-wall and overlap fermions 
are thus made more concrete. 

In order to develop this physical picture, we probe our 
quenched ensemble with the two-flavor fermion action 



ip(D - (A + jT 3 mi)7 5 )V' 
■ip(H w - (A + iT 3 mi))^', 



(1) 



where ip' = 75 tp, and Tk are Pauli matrices acting in flavor 
space. H\y = -D75 is hermitian. Explicitly, 



D = 



W + m -C 
W + m 



(2) 



where C xy = | Y?n=i {^x+p,,yU xll - Sx-ft,vU$ M ) and 
W xy = A5 xy ~ 5X^=1 ($x+p,,yU XIJ , + S x -p, tV U^). Each 
entry is a 2 x 2 matrix, with <r M = («?, i), and three 
Pauli spin matrices. The link variables U xfl £ SU (3) con- 
stitute the random field in which the fermions move. The 
parameter A will allow us to study the spectral density 
p{\) of Hw via a condensate, mi is a "twisted mass" 
which breaks isospin 0, 0] , and acts as an external mag- 
netic field for the condensate of interest. Neither A nor 
mi appears in the Boltzmann weight of the ensemble. 
For any A and mi one derives the Ward identity 

E d t < W *~(f)>+2roi (n + (x) 7r-(y)) = S xy fa(y)) . 

(3) 

Here d*f(x) = f(x) — f(x — jl) and J^(x) is the flavor- 
changing vector current jlfij . and tt^ = itpj 5 T±ip and 
7T3 = iip"f5T3tp , with t± = (ti ± iTi)j2. Introducing the 
Green function G = (Hw — A — imi)^ 1 one has 



<7r 3 ) = (2/V4)Imtr (G) 



(4) 
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where V4 is the four- volume. This implies a generalized 
Banks-Casher relation 

Urn (7T 3 ) = 2ttp(A). (5) 

mi— >0 

Thus the spectral density p(A) is an order parameter for 
flavor symmetry breaking in the quenched theory with 
fermion action Q). The easiest way to calculate p(A) is 
in fact through Eqs. (0J and JSJ. 

The two-point function T(x, y\ A) = (7r + (a;) ir~(y)) rep- 
resents correlations of the eigenmode densities of Hw- 
This is readily seen from its spectral decomposition, 

r(x,y;\) = (j2*tix)*k(xh J—— 

\ kn 

where ty n is the eigenmode of Hw with eigenvalue X n . 
We calculate it at zero three-momentum, A) = 
(t^) -1 X^j?r(0, x, t, y; A), where V3 is the spatial vol- 
ume. The mobility edge A c is determined as the value of A 
where these correlations become long-ranged as mi — > 0, 
that is, when the large-i behavior of T(t; A) changes from 
exponential (0 < |A| < A c ) to power law (|A| > A c ) in 
this limit. Above A c one has extended modes, and the 
long-range density-density correlations play the role of 
Goldstone bosons for flavor symmetry breaking. Below 
A c there are no long-range correlations, and no mass- 
less pole in ( J+(x) 7r~(y)). How, then, can the Ward 
identity © be satisfied in the limit mi — > 0? The an- 
swer is that, when p(X) arises from exponentially local- 
ized modes, the quenched two-point function T(x, y; A) 
diverges as 1/mi in this limit Hi. In fact 0, 

+0(1). (7) 

As mi — > 0, the expectation value in Eq. is non-zero 
if and only if p(X) ^ 0. It thus provides a mechanism for 
saturating the Ward identity without Goldstone bosons. 

If the mobility edge is at A = 0, Goldstone bosons 
dominate the correlation function; hence we may take 
A c = to be the definition of the Aoki phase. 

We have determined the value of A c at several locations 
in the (/3, mo) plane. For each value of j3 we generated an 
ensemble of 120 quenched configurations using the stan- 
dard plaquette action. The four- volume was 16 , with 
periodic boundary conditions for all fields. Our measure- 
ments were mostly done at mo = — 1.5 [hopping param- 
eter k = (8 + 2too) _1 = 0.2], which is roughly the value 
used in domain-wall and overlap simulations. We mea- 
sured T(t; A) using random sources on time slices and 
t. We extracted a mass M = M(A, toi) from T(t; A) at 



TABLE I: Exponential falloff rate of the density-density cor- 
relator F(t; A) at P = 5.85, mo = —1.5, for mi — > 0. Errors 
are statistical only. /1 is the extrapolation of masses deter- 
mined at mi 7^ 0, while jl is extracted by extrapolation of 
miT(t; X) to mi = 0. 



A 


M 2 




M 2 


0.0 


2.42(6) 




1.9-2.3 


0.1 


1.99(8) 




1.3-1.9 


0.2 


1.18(5) 




0.9-1.1 


0.3 


0.21(3) 




0.4-0.5 


0.4 


-0.04(2) 




0.12 


0.5 


-0.04(2) 




0.05 


0.6 


-0.05(2) 




0.01 


TABLE II: Mobility edge A c (/3,m ) 


and (when A c > 0) inter- 


polated spectral density p(A c ). Where no error is 


shown, the 


(statistical) 


error is less than one in 


the last digit. 







m 


A c 


P(Xc) 


00 


-1.5 


0.5 




6.0 


-1.5 


0.41 


0.14 


5.85 


-1.5 


0.32 


0.08 


5.7 


-1.5 


0.25 


0.07 


5.6 


-1.5 


0.14(2) 


0.05 


5.5 


-1.5 


0.0 




5.4 


-1.5 


0.0 




5.7 


-2.0 


0.21 


0.14 


5.7 


-2.4 


0.03(2) 


0.06 



toi values between 0.01 and 0.07, and extrapolated to 
toi = by fitting to M 2 = fi 2 (X) + a(X)m 1 17]. The 
results for (3 = 5.85 are shown in the second column of 
Table HJ One sees that fi 2 starts falling rapidly above 
A = 0.1. 

By definition /i 2 (A) drops to zero at the mobility edge 
A c . We determine A c by linear extrapolation from the 
last two points with positive /1 2 . Our results are compiled 
in Table ILT| Consider first the mo = —1.5 results. For 
reference, we include the free-theory limit {(3 = 00) 0], 
where A c coincides with the gap of the free Hw. At 
(3 — 6.0, A c is still close to its free-field value. The curve 
A c (/3) steepens before reaching zero somewhere between 
(3 = 5.6 and (3 = 5.5, where we enter the Aoki phase. 

Table ITT1 also shows results at two other mo values for 
(3 = 5.7. The mo = —2.4 result sug gests that one is near 
the boundary of the Aoki phase [l8| ■ We find only a small 
change between m = —2.0 and mo = —1.5, consistent 
with the finding of Ref. y that the spectral properties of 
Hw vary slowly over this range. This is why we have 
explored mainly the /3-dependence at fixed mo = —1.5. 

In order to account for the absence of a mass- 
less pole in Eq. Q, the 1/mi-divergence in T(x,y;X) 
must persist for a range of momenta p, and its co- 
efficient should depend smoothly on p. To confirm 
this, we calculated the Fourier transform F(w ra ;A) = 
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TABLE III: Spectral properties for (3 = 5.85, m = -1.5. 
The mobility edge is at A c = 0.32, marked by the horizontal 
line in the table. 



A 


P(A) 


R 


Is 


h 


0.0 


0.0011(1) 


17. 


3.1(4) 


0.64(1) 


0.1 


0.0019(1) 


15. 


3.5(5) 


0.71(1) 


0.2 


0.0088(4) 


10. 


3.7(4) 


0.92(2) 


0.3 


0.056(1) 


6.5 


4.3(3) 


2.2(2) 


0.4 


0.168(7) 


4.9 


4.7(8) 




0.5 


0.27(1) 


4.4 


11.(3) 




0.6 


0.39(2) 


4.0 


7.(1) 





J2t cos(ui n t)T(t; A), where u> n — 2irn/l6, and extrapo- 
lated m 1 f(tj„;A) linearly to mi = 01- Results are 
shown in Fig.UJfor A = 0.0 at (/?, m ) = (5.7, -1.5). The 
w-dependence of the 1/mi-divergence is indeed smooth. 
For comparison, we repeated the calculation for A = 0.5, 
which is above A c . The extrapolation of mir(0; A) to 
mi = is straightforward, as it must be since this gives 
p(X) according to Eqs. Q and J5J). Doing the same with 
r(w n ^ 0; A), however, leads to a huge x 2 - This confirms 
the qualitative difference between |A| < A c and |A| > A c . 

We present in Table II I II some quantities that further 
illustrate properties of the localized modes, for the same 
(/3,mo) as in Table [I] Using them, we can address the 
question of whether below the mobility edge p(X) arises 
from well-separated, exponentially localized eigenmodes. 

We define a generalized "participation ratio" Vnjpi a 
single eigenmode via V~ l = Et(J2s\^n(x,t)\ 2 ) 2 Q. If 
\^,i\ 2 has support mainly on a four- volume l^, then V n ~ 
l n . A spectral decomposition of the quantity V~ 1 (X) = 
lim mi ^o TO ir(i = 0; A) shows that it is equal to p(X) 
times an average of V~ x over eigenmodes with eigenvalue 
A„ = A. Thus, if we dehne the "support length" l s = 
p(X)P(X), we see that l/l s is an average of l/l n - The 
fourth column of Table Hill gives this l s , a measure of the 
linear size of the support of the eigenmodes. 

We may now compare l s to the distance between eigen- 
modes. We have a measure of the latter from the values 
of mi that wc used in measuring p(\). Spectral sums as 

0.01p ' 1 i 1 i 1 ' 1— 

0.008- J 
0.006- I 
0.004- 



0.002- 




FIG. 1: Coefficient of the 1/mi-divergence in T(a; n ;A) for 
P = 5.7, m = -1.5, A = 0.0. 



in Eq. J7J show that mi is the resolution with which we 
detect eigenmodes near A. For mi ~ 0.01, the number of 
modes we detect for a typical gauge configuration is thus 
N ~ 0.01%>(A), and so R(X) = (0.01 p(A))- 1 / 4 is a mea- 
sure of the average distance between modes. If l s -C R 
the modes are isolated, and correlation functions reflect 
properties of individual localized modes, with no inter- 
ference. Table IIIII shows this to be the case well below 
the mobility edge. 

The density of an exponentially localized mode has the 
asymptotic behavior 

|*(x)| 2 ~cexp(-\x-x \/li), (8) 

which defines k , the localization length. When the modes 
are isolated, the decay rate (extrapolated to mi = 0) of 
r(t; A) reflects the localization length of the individual 
modes. We thus define an average localization length 
through li = 1/fi (Table ITTT1 last column; compare Ta- 
ble[IJ 2nd column). Well below the mobility edge h turns 
out to be much smaller than l s . In fact, for j3 > 5.85, 
we find that k < 1 and l s > 3. This is good news for 
domain- wall and overlap simulations (see below). 

Equation (jSJ) represents an exponential envelope that 
we expect to multiply oscillations in | | 2 . These fluctu- 
ations, often large, survive the extrapolation of the cor- 
relation function T(t; A) itself to mi = 0. As a result, if 
we extract a mass p from lim mi ^o miT(t; A), the result 
varies with the details of the fit. We present rough values 
of p, in the last column of Table HJ In the absence of a 
model for the fluctuations, p is only a qualitative mea- 
sure, and hence we do not quote an error, p follows the 
trends shown by p. 

Taking all our results together, we have compelling evi- 
dence for exponential localization well below the mobility 
edge. Here the modes are isolated in the sense of l s <C R 
(for mi ~ 0.01), and k = 1/p provides an accurate esti- 
mate of the average localization length. For A > A c inter- 
ference effects destroy this connection. Upon comparing 
our data for all values of (/?, mo) shown in Table ITT1 we 
find that we can characterize the mobility edge itself as 
follows: (1) p(X) w 0.05-0.15 near A c ; (2) h> 1 signals 
the proximity of A c ; (3) / s ~ R ~ 5 at A c [llj . 

Finally, we revisit the implications of our results for 
domain- wall and overlap fermions. These two closely re- 
lated descendants of Wilson fermions employ a super- 
critical Wilson operator as a key element in their con- 
struction. Both are expected to be local, and they both 
have a (modified) chiral symmetry ^1 & t non-zero lat- 
tice spacing. The question is to what extent these ex- 
pectations are fulfilled in actual lattice QCD simulations 
employing these fermions. In Ref. we argued that lo- 
cality and chirality will coexist in these formulations if 
and only if the following holds: On a given ensemble of 
configurations, the mobility edge of the underlying Wil- 
son operator must be well above zero. It does not matter 
whether the ensemble is quenched or is generated with 
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dynamical domain- wall or overlap fermions. We will not 
repeat the whole argument Q here but rather focus on 
assessing the implications of our numerical results. 

Domain-wall fermions employ an auxiliary, discrete, 
and (in practice) finite fifth dimension with spacing 05 
and N s sites. Finiteness of the fifth dimension ensures 
locality but leads to "residual" violations of chiral sym- 
metry. A common measure of these violations, denoted 
as m ros , may be thought of as an additive correction to 
the quark mass. It is determined from a ratio of pseu- 
doscalar correlation functions at zero spatial momentum 
and time separation t. In Ref. it was argued that 



ci exp 



+ c 2 exp 



(ji(O)" 



(9) 

where the two terms arise from extended and localized 
modes, respectively, m v is the pion mass in the simula- 
tion. The mobility edge A c and the localization length 
k(X = 0) refer to a "hamiltonian" H obtained from the 
transfer matrix in the fifth dimension, which depends on 
05. Thus m ros reflects the spectral properties of H. We 
recover Hw from H in the limit 05 — > 0. Moreover, it 
can be proved that H^f = if and only if Hw^ — 0, i.e. 
the zero modes of H remain unchanged as 05 is varied. 
This implies that k(0) = k(0). For quenched simulations 
at ft > 5.85 we thus find that [;(0) = k(0) ~ 0.6. The 
last term in Eq. therefore vanishes rapidly with t. 
This agrees with previous findings ^(J that m rcs is fairly 
i-independent once t is large enough. Similarly, A c = 
if and only if A c = 0. Since good chiral symmetry re- 
quires TO rC s to be small, it follows that simulations must 
be performed well outside the Aoki phase of the under- 
lying Wilson operator. 

For overlap fermions, chiral symmetry is guaranteed, 
but not locality. Deteriorating locality may distort phys- 
ical predictions in an uncontrolled way. Indeed, for large 
separations one expects D ov (x, y) ~ c ov exp(— \x— y\/l ov ). 
The exponential tail of the overlap may be represented 
as an unphysical field of mass l/l ov that mixes with the 
physical quarks with an amplitude controlled by c ov - The 
range of the overlap operator, far from being merely a nu- 
merical nuisance, is thus a key indicator of the validity 
of a simulation. 

If an admissibility condition is imposed, it can be 
proved that the range Z Q v of the overlap operator is 0(1) 
in lattice units |13|. For realistic ensembles, / ov depends 
on the spectral properties of H\y, and good locality again 
requires keeping away from the Aoki phase. One an- 
ticipates that Z ov is on the order of either A^ 1 or k(0), 
whichever is larger; if it is the latter, c ov should be re- 
lated to p(0) 0. The spectral properties studied in this 
Letter are thus of central importance for understanding 
the locality properties of the overlap operator. 

Our computer code is based on the public lattice gauge 
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The factor mi/[(A n — A) 2 + mf] in Eq. Q justifies the 
linear fit. The contribution of modes with |A n — A| ~ mi 
to miP(t; A) is roughly constant, while that of "bulk" 
modes with |A n — A| S> mi vanishes linearly with mi. 



